A first principles study on FeAs single layers 
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FeAs~ single layer is tested as a simple model for LaFeAsO and BaFe2 As2 based on first-principles 
calculations using generalized gradient approximation (GGA) and GGA-|-i7. The calculated single- 
layer geometric and electronic structures are inconsistent with that of bulk materials. The bulk 
coUinear antiferromagnetic ground state is failed to be obtained in the FeAs~ single layer. The 
monotonous behavior of the Fe-As distance in z direction upon electron or hole doping is also in 
contrast with bulk materials. Our results indicate that, in LaFeAsO and BaFe2As2, interactions 
between FeAs layer and other layers beyond simple charge doping are important, and a single FeAs 
layer may not represent a good model for Fe based superconducting materials. 
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I. INTRODUCTION 

The recently discovered high temperature supercon- 
ductivity in LaFeAs[Oi_a;Fa,}i has attracted a lot of in- 
terests and triggered the research for other iron-based 
superconductors. Up to now, most of the Fe-based su- 
perconductors are based on either i?FeAsO (i?=La, Ce, 
Sm, Nd, Pr and Gd) or ylFe2As2 (A=Ba, Sr) structures. 
In the former 1111 series, replacing La with other rare- 
earth atom increases the transition temperature (Tc) up 
to 55 K for SmOi_a;Fa;FeAs:2-^ For the latter 122 series, 
with appropriate alkali metal (K and Cs) doping, Tc can 
be raised up to 37 Ki^i^ Very recently, superconductivity 
has also been reported for As-free material FeSei_5 and 
Fe(Tea;Sei_a;)i_5 with T^ around 27 K under pressurei^ 
Moreover, replacing RO layer in 1111 materials with Li 
or Na also leads to Tc of 18 K and 9 K respectively.*'^ 

The Fe-based superconductors have a quasi two- 
dimensional tetrahedral structure, where FeAs layers are 
separated by RO (i?=La, Ce, Sm, Nd, Pr and Gd), 
A {A=Ba, Sr), or Li (Na) layers. Except LiFeAs and 
NaFeAs, both parent compounds of 1111 and 122 super- 
conductors are metallic but not superconducting. They 
undergo a phase transition from tetragonal to orthorhom- 
bic with the decrease of temperature, which accom- 
panies with a new coUinear antiferromagnetic (AFM) 
order, known the SDW phase ;^°i^^'^^i^^ Upon doping 
the SDW is suppressed, and superconductivity appears. 
First principles calculations for both 1111 and 122 
materials have been reported, using either local spin 
density (LSDA) or generalized gradient approximations 
{GGA)^±^^^^^^^2J^ The coUinear AFM ground state in 
the parent compounds has been confirmed by theory. 

Although the mechanism of superconductivity in these 
Fe-based materials is still unknown, it is clear the essen- 
tial physics lies in the common FeAs layer. It serves as 
the conducting layer, and the interplay between mag- 



netism and superconductivity happens in this layer. 
Therefore, it is important to investigate its structure 
change as well as the evolution of electronic properties 
upon doping. 

At low temperature (25 K), in F doped LaFeAsO , 
with F doping, the Fe-As bond length changes less than 
0.1%, while the La-As distance reduces by ~1.5% and 
the La-0 distance increases by ^0.8%.~ These results 
demonstrate that the structure of FeAs layer changes 
slightly upon doping, in contrast with the significantly 
modified LaO layer. Besides, the layered structure of Fe- 
based superconductors is very similar to that of cuprates 
and Naa;Co02 superconductors. A single C0O2 layer has 
been successfully used as a model system to investigate 
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FIG. 1: (Color online) The crystal structure of the y/2 x 
\/2 X 1 single FeAs layer. The red and yellow balls are Fe 
and As species respectively, (a) is the side view of single 
FeAs layer, (b) and (c) are the top view of FeAs layer, where 
the black arrows signify the FM (b) and coUinear AFM (c) 
configurations on Fe atoms. 
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the doping effects on Naa;Co02'2S-'^ Thus, one question 
comes out, can we use a similar model of single FeAs layer 
to study the doping effects on Fe-based superconductors? 

In this article, we calculate the geometric, electronic, 
and magnetic properties of a singe FeAs~ layer in the 
framework of density functional theory (DFT). We fail 
to obtain the coUinear AFM phase with both optimized 
structure and experimental structure. The behavior of 
dAs (Fig- Ht') respecting to the doping level also differs 
from that of 1111 and 211 materials. 



II. MODEL AND METHOD 

As shown in FiglU FeAs layers are formed by edge- 
shared FeAs4 tetrahedras with Fe ions sandwiched be- 
tween two As sheets. In the undoped parent materials for 
both the 1111 and 122 series, the FeAs layer is negatively 
charged with one electron, which means undoped FeAs 
layer corresponds to FeAs~. Starting from this point, 
hole or electron doping is realized by simply adding or 
removing electrons from FeAs layer, with uniform com- 
pensated charge background. Doping level x is defined 
for charged system FeAs~'^^+^). In this work, doping level 
X ranging from —1.00 to -f 1.00 is investigated. 

Most of the theoretical works on the electronic and 
magnetic properties of Fe-based superconductors in the 
literature are based on the high-temperature tetrahedral 
structure, with a=6.— li^ii^ii^i^^ One attempt to obtain 
the low-temperature orthorhombic lattice constants used 
an alternative strategy with the magnetic moments fixed 
to experimental value instead of optimizing on ground 
state potential energy surface.— For single layer in our 
case, if we directly scan the orthorhombic structure pa- 
rameters, we obtain a and b around 4.64 A, which is sig- 
nificantly lower than the experimental values (by more 
than l.oA). Therefore, in the following calculations, we 
fix a and b to their experimental values 5.683 and 5.710, 
respectively*^^ In order to exclude interactions between 
neighboring layers, c is set to 14 A, corresponding to 
about 10 A's distance between two neighboring FeAs lay- 
ers. The positions of all atoms are allowed to relax until 
forces on each atoms are smaller than 0.01 eV/A. Due to 
the symmetry of the system, the only degree of freedom 
of the atoms is the z coordinate of As. As shown in Fig. 
[T^, z coordinate of As related to the Fe plane is marked 
as dAs- 

The electronic structure calculations are carried out 
using the Vienna ab initio simulation package.-^^ PBE 
functional is used,^^ The electron-ion interactions are de- 
scribed in the framework of the projected augment waves 
method and the frozen core approximation.— The energy 
cutoff is set to be 600 eV, the same as previously used 
for LaOFePi^ - For density of states (DOS) calculation, we 
used a 12x12x6 Monkhorst-Pack /c-point grid to sample 
the Brillouin zone, while for geometry optimization, a 
8x8x4 grid have been used. 

For magnetic property calculations, initial magnetic 




FIG. 2: Dependence of (a) dAs and (b) magnetic moment 
{fis/Fe) on Hubbard parameter U. The experimental value 
of dAs for LaFeAsO is marked by a dotted line. 



moments are set according to non-spin polarized (NM), 
ferromagnetic (FM), and coUinear antiferromagnetic 
(AFM) ordering. The FM and coUinear AFM configu- 
rations are illustrated in Fig[l] b and c. The latter is 
the experimentally observed ground state for parent com- 
pounds. 

In GGA-I-U calculations, we adopt a simplified model, 
where the on-site Coulomb repulsion U and the atomic- 
orbital intra-exchange energy J are simplified to one pa- 
rameter Ueff = U — J. For simplicity, we will call Ueff 
as U hereafter. 



III. RESULTS AND DISCUSSIONS 

First we performed structure optimization and elec- 
tronic structure calculations for FeAs~. From GGA re- 
sults, the total energy of FM state is about 0.002 eV 
lower than that of NM state, and the magnetic moment 
on Fe in FM state is about 0.18 fiB- The coUinear AFM 
state, which is reported to be the ground state of 1111 
and 122 materials, is not stable at all. 

The optimized dAs is 1.19 A, which is significantly 
lower than the experimental value for LaFeAsO (about 
1.31 A at 4 K).— The electronic structure as well as the 
magnetic moments on Fe in 1111 and 122 parent com- 
pounds are very sensitive to dAs so the discrepancy 
between the calculated magnetic structure and the ex- 
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FIG. 3: (Color online) Density of states of FeAs and the con- 
tributions of Fe d states and As p states to FeAs DOS. (a) 
PDOS of FeAs in LaFeAsO, (b) DOS of FeAs single layer with 
experimental values of Fe and As positions. The Fermi level 
is aligned to 0.00 eV. 

perimental one may cause by the heavily underestimated 

d-As- 

Considering the possible electron correlation in FeAs, 
GGA+U may be a useful strategy to correct (Ias- Dif- 
ferent U values are used to test it effects. In all DFT+U 
calculations, we get either an FM or a NM ground state. 
The magnetic moment increase with U for medium value 
of [/, and when U goes to relative large negative value, 
the magnetic moment on Fe will be quenched, which is 
consistent with the trend that negative U delocalize elec- 
tron. As shown in FiglS] (Ias decreases monotonously 
with U . To get the experimental (Ias value, an unphysi- 
cal U as low as -13.0 eV should be used. 

Since DFT-f U also fail to give the correct cIas , another 
thing we can try is to directly adopt the experimental dAs 
to see if we can get correct magnetic properties. The cal- 
culated total energy of FM state is about 0.0278 eV/cell 
lower than the NM state, with a magnetic moment of 
0.48 hb per Fe atom. The collinear AFM state does 
not exist again in single FeAs layer with experimental 
dAs- The resulting electronic density of states (DOS) of 
undoped FeAs layer is shown in FiglSl comparing to the 
FeAs partial DOS (PDOS) of LaFeAsO obtained with ex- 
perimental structure parameters at low temperature. ^^'^^ 
The main structures of the DOS and PDOS are very 
similar. In energy ranges from -2 eV to 2 eV near the 
Fermi level, Fe 3d states dominate. The peaks from -3 
to -2 eV are mixed Fe and As states. The density of 
states at Fermi level {N{Ef)) for single FeAs layer is sig- 
nificantly larger than that of FeAs PDOS of LaFeAsO 
(^^16. 13 eV~^ v.s. ~7.15 eV~^), which may result in 
different magnetic behavior. 
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FIG. 4: (Color online) The evolution of cIas with the doping 
level X. 

The doping effect for single FeAs layer is also stud- 
ied. Experimentally, the structure transition from tetra- 
hedral phase to orthorhombic phase is suppressed when 
sufficient electron or hole doping is applied, and the 
dAs in tetrahedral and orthorhombic phases in undoped 
LaFeAsO is almost the samci^ Therefore, the optimiza- 
tion for FeAs single layer with doping level from -1.0 
to 1.0 is carried out using the experimental tetrahedral 
phase lattice parameters (a=b=5.706 A). As shown in 
Fig. m the dAs changes monotonously with doping level, 
to be specific, increases with the level of electron doping, 
and decreases with the level of hole doping. In exper- 
iment, at 120 K for LaFeAsO, dAs do increase slightly 
from 1.319 to 1.323 A when 14% F is dopedi^i^o How- 
ever, in the case of hole doping, for Bai_a;K2,Fe2As2 at 
10 K, dAs increases from 1.344, to 1.351 A and 1.358 A 
when hole doping increases from zero to 10% and 20%, 
respectivelyiii»^ So the calculated trend of dAs using the 
single layer FeAs model is insufficient to describe the dop- 
ing effects on geometrical properties of Fe-based super- 
conductors. 



IV. CONCLUSION 

We have performed first-principles calculations on sin- 
gle layer FeAs. In this model where the inter layer inter- 
action is ignored, we find the structure of the FeAs layer 
in i?FeAsO and AFe2As2 can not be reproduced accu- 
rately in the framework of GGA and GGA+U. Besides, 
with both optimized and experimental lattice parame- 
ters, the collinear AFM ground state of i?FeAsO and 
AFe2As2 can not be obtained in FeAs single layer. 

In the simple single layer model, the inter-layer inter- 
actions between ReO (i?) layers and FeAs layers are ex- 
cluded. Our results suggest that this interactions may 
need to be considered to obtain correct geometry. This 
conclusion is important for choosing a proper theoretical 
model in future investigation of Fe based superconduc- 
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tors. 
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